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We study the effect of pure dephasing on the entanglement of a pair of two-level subsystems 
(qubits). We show that partial dephasing induced by a super-Ohmic reservoir, corresponding to 
well-established properties of confined charge states and phonons in semiconductors, may lead to 
complete disentanglement. We show also that the disentanglement effect increases with growing 
distance between the two subsystems. 



I. INTRODUCTION 

Entanglement 0, is one of the fundamental con- 
stituents of quantum theory. Correlations between the 
results of appropriately chosen measurements on entan- 
gled subsystems cannot be accounted for by any classical 
(realistic and local) theory Q , precluding the existence of 
a wide class of hypothetical more fundamental structures 
underlying the incompleteness of quantum description. 
Apart from its essential role in our understanding of the 
quantum world, entanglement is an important resource 
in quantum information processing J4| where it provides a 
quantum channel for teleportation |5| , superdense coding 
and distribution of cryptographic keys 0. 

In order to manifest genuinely quantum behavior re- 
sulting from entanglement a quantum system must main- 
tain phase relations between the components of its quan- 
tum superposition state, involving different states of dis- 
tinct subsystems. Keeping in mind that the subsystems 
may be separated by a macroscopic distance, one may ex- 
pect such a non-local superposition state to be extremely 
fragile to the dephasing effect of the environment. In fact, 
it has been shown that entanglement of a pair of two- 
level subs yste ms tends to decay faster than local coher- 
ence [§1 l9l Il0| . As expected, the decay of entanglement is 
stronger if the subsystems interact with different environ- 
ments (which might result from a large spatial separation 
between them): certain states that show robust entangle- 
ment under collective dephasing become disentangled by 
separate environments |9( . Remarkably, it was shown for 
different classes of systems QjJ, El ll3 that certain states 
may become separable (completely disentangled) within 
a final time under conditions that lead to usual, exponen- 
tial decay of local coherence. Since even partial entangle- 
ment of many copies of a bipartite quantum system may 
still be distilled to a smaller number of maximally entan- 
gled systems ^3|, many quantum information process- 
ing and communication tasks may be carried out using 
partly disentangled systems, as long as there is still some 
entanglement left. It is therefore essential to understand 
whether environmental influence leads to the appearance 
of separability in realistic models of dephasing. 

In this paper we study the decay of entanglement 
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between a pair of spatially localized two-level systems 
(qubits) under pure dephasing induced by a bosonic bath. 
In order to capture the physical aspects of disentangle- 
ment in real systems, we focus on the specific case of 
two excitons confined in one or two semiconductor quan- 
tum dots and coupled to phonons. In this system, the 
properties of the coupling and the resulting kinetics of 
dephasing are well understood. In particular, the ex- 
perimentally observed evolution of the system may be 
accounted for by pure dephasing within an independent 
boson model U4[. In contrast to the properties of the 
Ohmic model [la ] , the carrier-phonon dynamics in a real 
system leads only to partial dephasing 16]: non-diagonal 
elements of the density matrix do not vanish exponen- 
tially but rather decay only to a certain finite value (apart 
from other processes on much longer time scales). This 
behavior is characteristic of super-Ohmic spectral den- 
sities resulting from the actual carrier-phonon 
couplings and phonon density of states The effect 
of pure dephasing on a system of two excitons in a single 
quantum dot was studied in a recent work poj in the con- 
text of ultrafast optical excitation. Here we describe the 
dynamics of two spatially separated systems, focusing on 
the resulting loss of entanglement. 

The main result of this paper is that a physical mech- 
anism, based on a microscopic model of interactions be- 
tween charges and phonons in solids and quantitatively 
confirmed by experiments on real quantum dot systems, 
may lead to the complete decay of the entanglement of a 
pair of two-level systems. This is in striking contrast with 
the partial asymptotic decay of local coherence under 
the same environmental dephasing. Moreover, using the 
physical coupling constants for spatially localized states 
we are able to describe the effect of spatial separation on 
the evolution of entanglement and to physically account 
for the crossover from a collective to an individual reser- 
voir regime. These results are of importance not only for 
the general understanding of the properties of entangle- 
ment of open systems but also for practical tasks related 
to coherent control of charge states in multi-partite semi- 
conductor systems (like carrier states in coupled quan- 
tum dots). This, in turn, may affect the feasibility of 
solid-state based quantum computing schemes [^J |2^| , 
including quantum error correction based on collective 
encoding of logical qubits (concatenation) [23| . 

The paper is organized as follows. In Sec. [H] we define 
the model of the excitonic two-qubit system and find its 
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evolution. Next, in Sec. IIIII we discuss the evolution of 
entanglement between the two qubits. Sec. II VI concludes 
the paper with final remarks. 



II. THE MODEL AND EVOLUTION OF THE 
TWO-QUBIT SYSTEM 

In this Section we describe the evolution of the two ex- 
citonic qubits under the dephasing effect of phonons. The 
result will be used in the next Section for the calculation 
of entanglement between the two qubits. 

The system is described by the Hamiltonian 

H = ei(|l)(l|®I)+e 2 (l[® |l)(l|) + Ae(|l)(l|® |1)(1|) 

+(ii)<ii®i)£/£ i) (&i +&_„) a) 

k 

+(I ® |1)<1|) J2 fk\b{ + b- k ) + J2 ^b%, 

k k 

where the two states of each subsystem are denoted by 
|0) and |1), I is the unit operator, €1,2 are the transition 
energies in the two subsystems, Ae is an energy shift 
due to the interaction between the subsystems, f^'^ are 
system-reservoir coupling constants, bk,b\. are bosonic 
operators of the reservoir modes, and u>k are the corre- 
sponding energies (we put K = 1). The explicit tensor 
notation refers to the two subsystems but is suppressed 
for the reservoir components. 

Exciton wave functions will be modelled by anisotropic 
Gaussians with the extension Z / h in the xy plane for the 
electron/hole and l z along z for both particles. Then, the 
coupling constants for the deformation potential coupling 
between confined charges and longitudinal phonon modes 
have the form /£ L ' 2 ' ) = fke ±lkzd / 2 , where d is the distance 
between the subsystems (d — corresponds to a biexci- 
ton, i.e., two excitons with different spins in a single QD) 
and / fc = /( c) -/f } , with 
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Here V is the normalization volume of the bosonic reser- 
voir, fcj_/ z are momentum components in the xy plane 
and along the z axis, a e /h are deformation potential con- 
stants for electrons/holes, c is the speed of longitudi- 
nal sound, and g is the crystal density. In our calcu- 
lations we put <7 = 8 eV, cr^ = — 1 eV, c = 5.6 nm/ps, 
g = 5600 kg/m 3 (corresponding to GaAs), and l c = 4.4 
nm, Zh = 3.6 nm, l z = 1 nm. 

Although the above coupling constants correspond to 
a specific system, their form reflects rather general phys- 
ical conditions. First, for any translationally invariant 
reservoir, k may be interpreted as momentum. Second, 
the system is assumed to be localized in space and to in- 
teract with the reservoir via a spatially local interaction 
(this is natural for solid state systems and corresponds to 



the dipole approximation of electrodynamics). Then, if 
the system size is I the uncertainty of momentum is ~ l/l 
and momentum conservation allows coupling to bosonic 
modes only within this range of A;, as manifested by the 
Gaussian cutoff in our formula. The position-dependent 
phase is also a general feature, while the analytical be- 
havior at k — > is specific to the coupling. 

The Hamiltonian Q is diagonalized by the transfor- 
mation W = J2i=i where we use the standard 
product basis |0) = |0)|0), jl) = |0)|1), etc., and W{ are 
Weyl shift operators 



Wi = exp[ 



^9 k 
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unitary transformation one gets 



g k . Upon this 



where 



H = Wi?W f = H L + Hi + H I{ 



F L = £ 1 (|1X1|®I) + £ 2 (I®|1)(1|) 



describes the independent (local) evolution of the sub- 
systems, 

H I = AE(\1)(1\®\1)(1\) 
describes their interaction, and 



k 



is the reservoir Hamiltonian. The energies here are 



(i)|2 



and 



Ei = e t - 2^^fc|3fc 

k 
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The evolution operator generated by the Hamiltonian 
||TJ may now be written as U t = WWtUt, where U t = 
exp{-iHt) and W t = U t WU} . Since H is diagonal, the 
evolution described by Ut is trivial. The final formulas 
are further simplified by performing the local unitary ro- 
tation U L = cxp(iH L t) (note that [U h ,W t ] = 0). Then, 
for the reduced density matrix of the two-qubit system 
one finds p(t) = U~lp(t)U-L, with 

p(t) = Tr R [W+Wte-^ 1 *^ ® p T )e iHlt Wlw\ , (2) 

where we assumed that the system is initially in a product 
state with the reservoir in the thermal equilibrium state 
Pt- Since local unitary transformations do not change 
the amount of entanglement, we may use the density ma- 
trix p{t) instead of p(t) in the calculations. 
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The elements of the density matrix p~(t) are found using 
Eq. (J2J), the definition of the operator W, and rules for 
multiplying and averaging Weyl operators 0, 0] • The 
result may be written in the operator sum representation 



with the unitary operator 



(3) 



U = 



+ exp 
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and the set of Kraus operators 

K = diag[o(t),6(t),6(t),a(t)], 



Ki = diag[(a 2 (i) - l)Va 2 (t) + 1,0,0,0], 
K 2 = v/l- a 2 (t) diag[-a 2 (t), 0,0,1], 



K 3 = 
K 4 = 

where 

a(t) = exp 
b(t) — exp 



diag[0, 0, (b 2 (t) - 1) V& 2 (*) + 1, 0] , 
v/l-6 2 (t)diag[0,l,-fo 2 (;),0], 



k z d 



bfcl 2 cos 2 -^-(coswfet- l)(2n fc + 1| 
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and = fk/^k- The set of operators -ftf^ corresponds to 
a kind of two-qubit phase damping channel. This channel 
is generated by the same physical pure dephasing process 
that would lead to the usual phase-damping channel for a 
single qubit [2^1 . The operators K\p may be interpreted 
as an effect of a "charge detector" which is sensitive only 
to the total number of excitons in the system. On the 
other hand, ^3,4 reflect the action of a "discriminator" 
which detects in which of the two dots the exciton is 
present. Since the two subsystems are identical they may 
be distinguished by the reservoir only because of their 
different position. Therefore, the latter contribution to 
dephasing is inefficient for d — 0. For a system of two 
excitons confined in a single dot (d = 0) the present result 
reduces to that found in Ref. |20| . 

The dephasing action of the reservoir develops in the 
course of the joint carrier-phonon evolution. Initially, all 
the Kraus operators except for Kq are null. For long 
times, the factors cosujkt in Eqs. i)4alb) become quickly 
oscillating functions of k and their contribution averages 
to 0. Consequently, the operators reach a certain 
asymptotic form. As a result, the non-diagonal elements 
of the density matrix decrease form their initial value 
to a certain asymptotic value depending on material pa- 
rameters, system geometry and temperature (see Fig.^l. 
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FIG. 1: The evolution of a two-qubit state under phonon- 
induced pure dephasing. The plot shows the relative reduc- 
tion of the non-diagonal elements of the reduced density ma- 
trix, nj = \pij(t)/pij(0)\ for T = 40 K, d = (left), and d = 6 
nm (right). For identical qubits one has roi = ro2 = rig — T23- 



For the actual carrier-phonon couplings in a semiconduc- 
tor system this asymptotic value is nonzero, similarly as 
in the single-qubit dephasing and t ne phase informa- 
tion is only partly erased. The non-monotonous evolution 
of coherence is due to the fact that the process is non- 
Markovian and some coherence may be regained in the 
course of system-reservoir interaction before the process 
is completed. It is also clear from Fig. ^ that the coher- 
ence between the states |1) = |0)|1) and |2) = |1)|0) may 
be broken only if the reservoir can distinguish between 
the two systems (i.e., for d 7^ 0). 

In order to interpret the time dependence of the ma- 
trix elements for d 7^ 0, shown in Fig. ^ (right), let us 
recall that the physical mechanism of dephasing is the 
emission of phonons from the excitonic qubit to the bulk 
of the crystal pCil [271 • These phonons form a spherical 
wave packet around the quantum dot, extending with 
the speed of sound and thus carrying the information 
from the qubit to the outside world [24| . For short times 
t < d/c this reservoir perturbation pertaining to each 
qubit stays localized around this qubit (the phonon wave 
packets do not overlap) so that detecting the presence of 
a charge is as efficient as deciding in which dot it is local- 
ized. As a result, the coherences between the ground and 
two-exciton state (matrix element P03) and between the 
two single-exciton states (P12) are affected to the same 
extent. For longer times, the wave packets originating 
from the two qubits intersect which, on one hand, leads 
to positive interference and to increased efficiency of over- 
all charge detection. This is reflected in the right panel 
of Fig. n°y a step- wise drop of po3 at t ~ d/c. On the 
other hand, however, this partial overlap makes it harder 
to identify the origin of the perturbation, hence the "dis- 
criminator" measurement becomes less efficient. There- 
fore, at t > d/c the inter-qubit coherence p\2 becomes 
less affected. 



III. DISENTANGLEMENT BY PURE 
DEPHASING 

For a quantitative description of the decay of entangle- 
ment, a measure of entanglement that may be calculated 
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from the system state is needed. For pure states, the von 
Neumann entropy of one subsystem [2g may be used but 
for mixed states there is no unique entanglement measure 
[2^,l30| . In order to quantitatively describe entanglement 
of systems in mixed states one defines various quantities, 
which may give essentially different entanglement mea- 
sures. Nonetheless, any such quantity must satisfy cer- 
tain natural requirements in order to be useful as an en- 
tanglement measure: it must vanish for separable states 
and be non-zero for entangled ones; it must reduce to 
the von Neumann entropy for a pure state; it must be 
non-increasing under local operations (it is impossible to 
increase entanglement by manipulating only one subsys- 
tem). One choice is to use the entanglement of formation 
(EOF), defined as the ensemble average of the von Neu- 
mann entro py m inimized over all ensemble preparations 
of the state |29l |3(j i an d thus being a natural generaliza- 
tion of the von Neumann entropy to mixed states. Such a 
measure may be interpreted as the asymptotic number of 
pure singlets necessary to prepare the state by local op- 
erations and classical communication. A practical char- 
acterization for mixed state entanglement is available for 
small systems but an explicit formula for calcu- 

lating an entanglement measure is known only for the 
EOF of a pair of two-level systems [^3, 0| ■ 

We will perform the calculations for two initial pure 
states 



l4 1} > 



|00) + |01) + |10) - |11) 



(5a) 
(5b) 



It is assumed that these states are prepared within a 
time shorter than reservoir response times. Both initial 
states l|5afl and (|5b|) are maximally entangled, i.e. their 
entanglement of formation is equal to 1. 

Using Eq. J3J), we calculate the EOF of the system 
state at a time t from the Wootters formula 33, 341 



E[p{t)] = E[p(t)} = -x+ log 2 x + - x^ log 2 x. 



where 



l±^l-C2{p(t)} 



and C[p{t)\ is the concurrence, given by 

C[p(t)} = max(0, A - Ai - A 2 - A 3 ), 



(6) 



where Ai are the eigenvalues of the matrix p(t){o~ y ® 
CJ y)p{t){ (J y ® a y) m decreasing order. 

The evolution of the EOF of the qubit pair is shown 
in Fig. |2 In the absence of energy shift AE (solid lines) , 
entanglement decays on a time scale of a few picosec- 
onds. At low temperatures or for overlapping systems, 
this process resembles the decay of coherences shown in 
Fig. n However, for a sufficiently large separation be- 
tween the systems and at sufficiently high temperatures 




FIG. 2: Evolution of entanglement of the two-qubit system 
at various temperatures with d = 6 nm (a,b) and for various 
distances d at T = 100 K (c,d). The left panels (a,c) show 
the result for the initial state 15aH and the right ones (b,d) 
for the singlet state (I5bt . The inset in (c) shows the lowest 
temperature at which the decay becomes complete for a given 
distance. Solid lines correspond to AE = and the dashed 
line to AE/h = 6 ps" 1 (at T = 40 K). 



the initially maximal entanglement present in the state 
(|5a|l decays completely after a finite time even though 
the environment-induced dephasing is always only par- 
tial (see Fig. J2Ji,c). The temperature T c at which the 
system becomes separable is related to the distance as 
shown in the inset in Fig. On the other hand, for the 
other initial state [Eq. (|5b|) ]. the destruction of entangle- 
ment is always only partial (Fig. Ej3,d). 

An important case is that of AE ^ 0. Such an energy 
shift (known as the biexcitonic shift in a semiconduc- 
tor system) leads to an entanglement-generating evolu- 
tion. This mechanism is used for performing nontriv- 
ial two-qubit gates (controlled-shift) in many proposals 
for semiconductor-based quantum information process- 
ing 0, . As can be seen in Fig. (dashed line) , in 
the presence of phonon-induced pure dephasing the cyclic 
evolution of entanglement is damped and the maximum 
achievable level of entanglement is reduced. Moreover, 
extended periods of time appear when the entanglement 
remains zero. 

The appearance of complete disentanglement for some 
initial states under sufficiently strong partial pure de- 
phasing may be understood with the help of Eq. JfjJ). If 
the completely dephased state (with a diagonal density 
matrix) has Ao — Ai — A 2 — A3 < then, by continuity, 
it will be surrounded by states with vanishing concur- 
rence, thus separable. In this case entanglement vanishes 
for sufficiently strongly dephased states, before the com- 
plete dephasing is reached. For a diagonal density matrix 
one finds A - Ai - A 2 - A 3 = -2min(pooP33,PiiP22), so 
that the above condition may only be satisfied if all four 
diagonal elements are nonzero. For instance, the totally 
mixed state (with a density matrix proportional to unity) 
is surrounded by a ball of separable states [3(| . Out of the 
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FIG. 3: The time at which the entanglement of the initial 
state (l-~>at decays completely as a function of temperature. 
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FIG. 4: The asymptotic value of the EOF as a function of 
the distance between the subsystems. Left and right panels 
correspond to the initial states I5all and 15bll . respectively. 

initial states considered here, the first one [Eq. Q5a[l] sat- 
isfies the above condition (it decays towards the totally 
mixed state) but the second one [Eq. (|5b|) ] does not. In 
general, the structure of entangled and separable states 
around the final state may be quite complicated. 

The time at which the entanglement of the state H5a|) 
vanishes completely depends on temperature and on the 
distance between the systems. As can be seen in Fig. |3| 
this time becomes finite only at a certain temperature 
(dashed line in the figure). Slightly above this criti- 
cal temperature, complete disentanglement takes place 
only for strongly separated systems. For higher temper- 
atures the disentanglement time for non-overlapping sys- 
tems depends very weakly on the distance. It should 
be stressed that the appearance of complete disentan- 
glement at increased temperatures is only related to 
stronger dephasing and, in principle, the state might be- 
come separable already at T = if the coupling were 
sufficiently strong. 

On the other hand, the dependence on the separation 
between the two subsystems reflects a more fundamental 
crossover, from the regime of a common reservoir to that 
of independent reservoirs 0- In our model, the non- 
diagonal element pi2 is unaffected for completely over- 
lapping systems (d = 0) (see Fig. [Q. This single non- 
vanishing element is enough to prevent complete dis- 



entanglement at any temperature and for any coupling 
strength. Only when the distance between the systems 
grows, dephasing is able to decrease this element to the 
extent sufficient for the total destruction of entangle- 
ment. This distance effect is shown in Fig. where we 
plot the asymptotic (long-time) value of the EOF. 

The effect of spatial separation between the systems is 
very clearly visible for the second initial state [Eq. I]5b[l]. 
which involves only this robust non-diagonal element (see 
Fig. 01 right). Here, the entanglement is absolutely stable 
if the systems overlap (see also Fig. |5Ji) but becomes 
much more fragile as soon as the separation between the 
systems is comparable to their size. This demonstrates 
that the distance between the subsystems is the physical 
parameter that governs the crossover between the two 
regimes of entanglement decay. 



IV. CONCLUSIONS 

We have studied the effect of pure dephasing on the 
entanglement of two-level subsystems. We have shown 
that partial dephasing may be sufficient to completely de- 
stroy entanglement for a class of initial states that may be 
easily characterized. Apart from its dependence on the 
initial state and temperature, the disentanglement effect 
shows essential dependence on the distance between the 
subsystems, manifesting a crossover between two regimes 
of reservoir-induced dephasing. Complete disentangle- 
ment appears only for spatially separated systems. 

Complete disentanglement due to partial pure dephas- 
ing typical for localized carrier states in semiconductor 
systems is not only of general interest but also of rele- 
vance to solid state implementations of quantum informa- 
tion processing. Moreover, understanding the role of the 
distance between subsystems in maintaining quantum 
correlations between them seems to be essential for re- 
alistic design of quantum error correction schemes based 
on collective encoding of logical qubits (concatenation) . 

One should note that, due to a purely non-Markovian 
character of the dephasing effect, the system evolution 
depends on the way in which the "initial" state has been 
prepared from the ground system state. Therefore, a 
reduction of the destructive effect may be expected if the 
preparation is done either slowly (adiabatically) or by 
shaped pulses [33, HU . 
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